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ABSTRACT 

The Mid Infra Red Instrument (MIRI) is one of the four instruments onboard the James Webb Space Telescope (JWST), 
providing imaging, coronagraphy and spectroscopy over the 5 — 28 /im band. Tj verify the optical performance of the 
instrument, extensive tests were performed at CEA on the flight model (FM) of the Mid-InfraRed IMager (MIRJM) at 
cryogenic temperatures and in the infrared. This paper reports on the point spread function (PSF) measurements at 5.6 /im, 
the shortest operating wavelength for imaging. At 5.6 /im, the PSF is not Nyquist-sampled, so we use am original technique 
that combines a microscanning measurement strategy with a deconvolution algorithm to obtain an over-resolved MIRIM 
PSF. The microscanning consists in a sub-pixel scan of a point source on the focal plane. A data inversion method is used 
to reconstruct PSF images that are over-resolved by a factor of 7 compared to the native resolution of MIRI. We show that 
the FWHM of the high-resolution PSFs were 5 — 10 % wider than that obtained with Zemax simulations. The main cause 
was identified as an out-of-specification tilt of the M4 mirror. After correction, two additional test campaigns were carried 
out, and we show that the shape of the PSF is conform to expectations. The FWHM of the PSFs are 0.18 — 0.20 arcsec, 
in agreement with simulations. 56.1 — 59.2% of the total encircled energy (normalized to a 5 arcsec radius) is contained 
within the first dark Airy ring, over the whole field of view. At longer wavelengths (7.7 — 25.5 /im), this percentage is 
57 — 68 %. MIRJM is thus compliant with the optical quality requirements. This characterization of the MIRJM PSF, as 
well as the deconvolution method presented here, are of particular importance, not only for the verification of the optical 
quality and the MIRI calibration, but also for scientific applications. 
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1. INTRODUCTION 

The JWST Mid Infra Red Instrument (MIRI, Wright et al. 2004,' Rieke et al. 2005^) is part of the scientific payload of 
the James Webb Space Telescope (JWST) (Fig.[T]i. Three models were built before the Flight Model (FM): the Structural 
Qualification Model (SQM), the verification model (VM), and the Engineering and Test Model (ETM). The different steps 
leading to the FM integration of the Mid-Infra Red IMager Optical Bench (MIRIM-OB), and the principal results associated 
to the SQM for vibration, and the ETM for optical performance, are presented in Amiaux et al. 2008.^ 

To verify the optical performance of the MIRIM FM, extensive tests were performed from Dec. 2008 to Dec. 2009 
at CEA, in the infrared and at cryogenic temperatures. An overview of the FM tests is presented in Ronayette et al. 
2010.13 In this paper we focus on the characterization of the point spread function (PSF) of MIRIM, and particularly at 
5.6 /im. This is the shortest operating wavelength, and therefore the most critical to assess the optical quality of the imager. 
We use a microscanning technique to measure and reconstruct the PSF at a higher resolution than the native one of the 
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Figure 1. Left: The Mid Infra Red Instrument (MIRI) as part of the scientific payload of the JWST. Middle: The MIRI Optical Bench 
Assembly (MIRI-OBA), with the imager in blue. Right: The MIRIM (Mid Infra Red IMager) mechanical layout, with the filter wheel 
assembly (FWA). MIRIM shall provide the following science functions: photometric Imaging between 5 and 27 fim, coronagraphy 
between 10 and 27 fim, low resolution spectroscopy between 5 and 10 fim. 



instrument. This allows us to characterize the MIRIM PSF very accurately and verify that the imager is compliant with 
PSF requirements. This study will be particularly helpful for the photometric calibration of MIRI, but also for scientific 
applications. 

We present three series of tests, resulting from the FMl (December 2008), FM2 (April 2009) and FM4 (December 
2009) cryogenic test campaigns at CEA. The field of view (henceforth FoV) and PSF measurements carried out during the 
FMl campaign showed a defect in the optical quality (vigneting issue and wider PSF than expected, see Sect. 16. 11 1. We then 
show the results of the FM2 and FM4 campaigns, after correction of this defect. This paper is organized as follows: Sect.|2] 
briefly presents the instrumental setup used to simulate the JWST beam and perform the cold tests. After a short description 
of the data reduction (Sect. |3]l, we derive the response curve of the detector and show how we correct the images for the 
non-linearity of the detector (Sect.|4|i. Then, Sect. |5] describes the microscanning test and the method used to reconstruct 
the MIRIM PSF at high resolution. Sect. |6]presents our results about the characteristics of the PSF. We finally present our 
conclusions in Sect. [T] 



The verification of the MIRIM optical quality has been done with an ambient temperature Telescope Simulator (TS), 
installed outside the helium-cooled cryostat that contains MIRIM-OB (Fig.|2]l. The TS simulates the optical beam delivered 
by the JWST. The optical conception of the TS is based on two off-axis parabolic miiTors (PMl and PM2) and a pupil 
mask in between. The pupil mask is the pupil of the telescope simulator (so-called the "STOP") and it reproduces the 
pupil of JWST. It is mounted on a motorised translation and rotation stages (4 axis: 3 translations and one rotation), which 
allows us to co-align the telescope simulator and the MIRIM pupils. The FoV is 80 mm x 80 mm for a MIRIM FoV of 
72 mm x 72 mm. An analysis of the image quality shows that the distortion amount of the TS is 0.33% for the useful field 
(Y. Longval, private communication). The aperture is i^/20 in the image plane. The exit pupil location is 3017.5 mm from 
the image plane for a diameter of 151.6 mm. The pupil is 75.4 mm in diameter To implement the MIRIM cryostat in a 
minimum space, folding mirrors are used. 

An IR point source with a shutter is mounted on a remotely controlled hexapod, that has 6 degrees of freedom (X, dx, 
Y , dy, Z, 9z), with an accuracy of 1 /im. The minimum diameter of the point source is 30 /^m, and the temperature of the 
source can be adjusted at 1 150 K or 2000 K. An extended black body source (400 K) can also be used for flat fields for 
instance. The lOK screen in the cryostat is fitted with a neutral density (ND) that reduces the flux due to the ambiant warm 
temperature (295 K) background. 



2. INSTRUMENTAL SETUP: THE MIRIM TEST BENCH AT CEA 
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Figure 2. The MIRIM test bench. The left panel shows the optical scheme of the Telescope Simulator (TS) that simulates the JWST 
beam and pupil to characterize MIRIM in image or coronographic mode. The warm TS is installed outside the cryostat (in blue on the 
right panel) that contains MIRIM. The TS allows us to use a point or extended source. 



For the microscanning, the shutter is used in chopping mode. Several frames are taken for each cycle (open— close shutter). 
For each cycle, the first frame is skipped because during this frame the shutter is closing. Then, the following 3 frames 
are skipped because of latency effects on the detector since the point source is very bright. For the FM2 test campaign, 
sub-arrays of 72 x 64 pixels were read out to speed up the data acquisition and reduction, with a detector integration time 
DIT ~ 0.86 s. Each exposure comprised 10 cycles, and the final image is obtained by taking a median from these 10 
cycles, corresponding to a total integration time of Tcxp = 132.3 s. For the FM4 campaign, the acquisition settings were 
different. 136 x 128 sub-aiTays were read out, with a detector integration time DIT — 0.30 s. The finale image is a median 
created from 50 cycles, corresponding to a total integration time of Tcxp = 183.3 s. The final images are also flat-fielded. 



The MIRIM detector (a Raytheon SB305 engineering model type) is a 1024 x 1024 pixels SiAs aiTay, loaned to CEA by 
JPL for the purpose of the test campaigns. The detector configuration is a Sensor Chip Assembly (SCA), i.e. an ultra-light 
mechanical and electrical housing. CEA has developed an appropriate housing (similar to the MIRI focal plane) to mount 
the IR detector onto the imager for the cold tests. Thus, qualitative analysis of the optical properties, including stray light 
performances, are possible. Note that the detector and the readout electronics used for these tests were not FM hardware. 

We tested the response of the SB 305 detector at CEA, by measuring the flux of an extended source with increasing 
integration times. The temperatures of the detector and the extended source are constant (5 K and 30 K, respectively), 
the integration time is the only parameter that changes. No filter has been used. Fig. [3] presents the curve of the detector 
response, i.e. the read value of the flux (in ADU) as a function of the "coiTected value". The coiTected value is the expected 
flux if the response were linear The 30 K blackbody source does not allow us to explore the linear part of the response 
curve, at low fluxes. With the smallest integration time (2.7 s for a full frame), we reach a signal of 26 360 ADU, which is 
already in the non-linear regime (Fig. [3]). To measure low fluxes in the linear regime, we used colder black body sources. 
An accurate sampling of the response curve at low ADU values shows that the response is linear for ADU values lower 
than w 16 000 (V. Moreau, private communication) and follows the law: 



3. DATA ACQUISITION SETTINGS AND REDUCTION 



4. LINEARITY CORRECTION FOR THE RESPONSE OF THE DETECTOR 



S'lin = a X tint + P , 



(1) 
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Figure 3. Left: response curve of the IVIIRI SB305 detector obtained at CEA. The baclc curve show the average read ADU values of the 
pixels in the center of the field as a function of the corrected values (see text for details). The red curve indicates an exponential fit (see 
§|4] Eq.[2l(. Right: impact of the linearity correction on the low-resolution PSFs. The plots show IVIIRIIVI PSFs profiles at 5.6 /im in two 
perpendicular directions. The pixel size is 25 /im. The black line represents data that have not been corrected for the detector response. 
The red dashed line shows the same PSF profile when the linearity correction is applied (see Sect.|4]l. The PSF shown here were taken 
during the FIVI2 campaign, in the center of the IVIIRIIVI field of view. 



where Sun is the signal in ADU in the linear part of the response, a = 3067 ADU s~^ is the slope of the linear response 
and /? = 18030 ADU the offset. Note that absolute photometric measurements were not possible because the quantum 
efficiency and conversion gain of the detector used for these tests were not characterized. 

Fig. |3] shows that the detector is non linear over a wide range of ADU values. The response curve can fitted with an 
exponential function: 

5,.ead = A-Se-«-/^, (2) 

where 5rcad is the read signal in ADU and 5cor the connected one. We found A — 56 418.5, B — 62 939.1 and C — 
35 990.5. 

To apply the linearity correction, we derive S'cor from S'rcad and correct all the individual images. To quantify the impact 
of the non-linear response curve on the PSF measurements, we provide two sets of test data: one that is not corrected for 
the detector response, and the other which is corrected for the non-linearity. We perform the PSF analysis on the two data 
sets and we compare the results. Fig. [3] compares the PSF profiles when the linearity correction is applied or not. The 
correction for the non-linear response of the detector improves the sharpness of the MIRJM PSFs because the deviation to 
the linearity is greater at high fluxes than the one at low signal. The FWHM of the corrected PSFs are ~ 5 % lower than 
the raw ones. Note that, although we do see a non-linear response in the FM detectors, the deviation to a linear behaviour 
is not as severe as in the detector used for these tests. 

5. MICROSCANNING TEST AND DATA INVERSION METHOD 

We describe here the aims of the microscanning tests and the deconvolution method we used to reconstruct over-resolved 
images of the PSF. The results are discussed in Sect. |6] 

5.1 Aims and method 

The FWHM of the MIRIM PSF at 5.6 ^m is 0.185", i.e. less than two pixels (1 pixel3= 0.11"), so the PSF is not Nyquist- 
sampled at this wavelength. The goal of the microscanning test is to obtain an over-resolved image of the PSF from 
observed multiple low-resolution (LR) images. Such a resolution enhancement approach has been an active research area 
(e.g. Park et al. 2003, El Molina et al. 2006,'^Idier 200^, and it is sometimes called over-resolution image reconstruction. 

The method consists in scanning a point source in a fine spatial resolution. The microscan pattern is given in Fig. |4] 
A 11 X 11 points grid was done, corresponding to a total scanning area of 1 MIRJM pixel on the detector plane, i.e., an 
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Figure 4. Left: Pattern of the microscan. All x 11 points scanning grid is performed, corresponding to a total scanning area of 1 pixel 
on the detector plane. The source is moved by steps of 7 /xm on the hexapod mount of the telescope simulator. This corresponds to 
individual displacements of 2.5 /im on the detector plane. Right: Sketch of the different steps to reconstruct the high-resolution image 
(inverse problem). 



area of 25^ /itn^. The source has been moved by steps of 7 /im on the hexapod mount. This corresponds to individual 
displacements of 2.5 /xm on the detector plane. The best focus was searched at the beginning of each test serie. Microscans 
were performed at 9 positions of the imager FoV. 

We summarize how we reconstruct the over-resolved images. Fig. |4] shows a sketch of the inverse problem. Let us 
denote x the high-resolution (HR) PSF. The aim of the deconvolution method is to deduce x from N low-resolution 
images yk (the observed data). Assuming that we start from the high-resolution image x, the forward problem can be 
expressed as 

Vk = SRTkX + Uk (3) 

where Tk, Ret S are three matrices associated with the operations of translation, convolution by the impulse-response of 
the detector (PSF), and down-sampling, respectively, k e {1, . . . , N}, where N is the total number of images taken at 
different positions on the detector n is the noise, whose level is assumed to be commensurate for all the images. Note that 
in Eq.|3] the only operation that is different from an image to another is the translation Tk- 

The first step to reconstruct HR images is to co-add the individual LR images on a fine grid, over-sampled by a factor /. 
To do so, the relative translations between the LR images have to be known with precision. Unfortunately, the mechanical 
control of the hexapod (where the point source is mounted) is not perfect, which results in mismatches between the expected 
position and the effective observed location of the source on the detector. Therefore we chose to estimate the shifts between 
images by cross-correlation. The three steps of the reconstruction are the following: 

1. we first estimate the translation between the low -resolution (LR) images by cross correlation (see Sect. l5.21 l. 

2. then, we co-add the LR images on a fine grid, assuming that there is no convolution by the PSF of the detector. 

3. finally, we solve the deconvolution problem with a Bayesian approach to estimate x (see Sect. l5.3l l. 

The deconvolution algorithm is a simplified version of a more complex algorithm developped by Rodet et al. 20051^ to 
perform over-resolution with spectral mapping data provided by the Spitzer telescope InfraRed Spectrometer (IRS). 

5.2 Estimate of the translations between images and co-addition 

We estimate the relative motions between the images by cross-correlation. To increase the precision of our estimate of sub- 
pixel translations between images, we perform a bilinear interpolation on an HR grid at a resolution equal to three times 
the over-sampling factor (3 x /). We define the first image as the reference image Jiof and we search, in an exhaustive 
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Figure 5. Expected vs. Estimated displacements of the source in the X and Y directions during the microscanning. A slight discrepancy 
is seen between the estimated moves on the image and the expected position of the source for the F-direction because, in practice, to 
control the absolute position of the hexapod is difficult. For the X-direction, the agreement is excellent. (Center of the FoV, FM2 data). 



manner, the shift that maximizes the cross-correlation factor, Cor, between the reference image and the current translated 
one, I. The cross-correlation is computed as: 

Cor(J,/,rf) = Yl E - W) (^-f-.j - (^^f)) (4) 

In Eq.m o-/^.^j. and aj are the standard deviations of the matrices /fef and J, respectively, (/ref) and (J) are the averages 
of /i-cf and I, respectively. 

Then, these motions are compared with the expected values measured with the hexapod. Fig. |5] shows an example of 
the comparison between the read position of the hexapod (expected position in blue) and the position estimated by cross- 
correlation (in red). The agreement between expected and estimated positions for the displacement in the X-direction is 
very good. However, a discrepancy is seen for the F-direction. The hexapod did not move at the expected positions in this 
direction. This issue is likely to be due to the control system of the hexapod support, which is not optimally designed to 
perform translations along the X and Y axis. 

The mismatch between estimated and expected relative motions between images introduces errors in the reconstruction 
method if one uses the expected positions. Therefore, the estimated motions by cross-correlation are used to reconstruct 
the high resolution image. Assuming that the translations Tk are known, we co-add low-resolution images yk onto an 
over-resolved grid. Thus, a HR image z is created, with a number of pixels times larger than the LR image, where 
/ is the over-sampling factor. However, this image is still bluiTed because of the response of the detector. The slight 
discrepancy between expected and estimated positions represents a lack of information, in particular in the F-direction. 
Some pixels of the HR image z are not present in the data yf^. 

Ideally, the over-sampling factor may be chosen up to a value of / = 10, but because of this lack of information, this 
choice introduces some slight artefacts in the reconstructed images. We chose a value of / = 7, which appeared to be a 
good trade-off between the minimization of artefacts and the gain in resolution. 

5.3 Deconvolution 

z is the convolution of the HR image x that we want to calculate by the impulse-response (PSF) of the detector: 



z = Rx + n 



(5) 



One cannot apply the classical methods of deconvolution because part of the data is not observed (due to uncertainties on 
the translations, the source is not exactly at its expected position on the fine grid image, see sect. 15.21 ). We take into account 
the truncation of the data by introducing a matrix A in Eq.|5] 

z — ARx + n = Hx + n (6) 

A is a mask with pixels value of 1 where the data exists, and otherwise. The operator H = AR returns the measurements 
z from the HR image x. 

To solve the deconvolution problem, a Bayesian formalism is used, and the lack of information in the data is supplied by 
the prior. In practice, the deconvolution consists of minimizing a least-square criterion. The estimate of the value x is given 
by the argument of the maximum of the posterior probability function, often called maximum a posteriori (MAP). Since 
the noise n is assumed to be Gaussian, independent, and uniform over the whole image, we use a Gaussian likelihood, and 
we obtain the following least-squares estimator: 

XMAP = arg min(||z - Aila;|p + fx\\Dx\\'^) (7) 

where fj, is the regularization parameter This equation defines the regularized least-square estimator, where D is the finite 
differences operator (Laplacian). The least-square criterion, Qls. to minimize in Eq.|2]is quadratic and can be re-written 
as: 

QLsix) = {z-HxY{z-Hx) + nx^D^Dx (8) 

The solution is explicit: 

XML = {H'H + fiD'D)-'H'z (9) 

The main issue with this explicit solution is that the size of the matrix H^H can be large and its inversion can be very 
time-consuming. Since the criterion to minimize is quadratic, it is convex and has a unique minimizor. We use a method 
called "gradient descent" or "steepest descent", which consists in starting from whatever value for x, and then converge 
to the solution in the direction of the steepest slope (opposite of the gradient's direction). At each step, this algorithm 
ensures that, after a sufficient number of iterations, we will approach the required minimum value of the criterion. Thus 
we calculate the gradient of the criterion to minimize: 

g{x) = ^9]^^ = 2 H\Hx -z)+ 2fiD'Dx (10) 
ox 

and we update the solution with the following equation: 

= _ aopt{H\Hx'' - z) + fiD'Dx'') , (11) 

where aopt is the optimal step to ensure the most efficient descent is the opposite direction of the gradient. This step is 
calculated by cancelling the gradient with respect to a. In our case, we use a variant of the gradient descent, so-called 
"conjugate gradient descent", which consists in having orthogonal directions at each descent steps. 



6. RESULTS: THE HIGH-RESOLUTION PSF OF THE MIRI IMAGER 

Fig. |6] compares the 5.6 pm PSF at the native resolution (LR, left image), with the over-resolved PSF obtained with the 
deconvolution method described above (middle image), and the simulated (Zemax) PSF at high-resolution (HR). The 
linearity correction is applied on the data and the over-sampling factor is set to / = 7. LR data correspond to 25 /im pixels 
on the detector, wheras HR correspond to a pixel size of 3.6 /im. One original MIRIM pixel equals to 7'^ pixels in the 
over-resolved image. 

The secondary Airy diffraction ring, which is barely visible on LR images, is clearly exhibited on the HR images. The 
shape of the secondary lobes is resolved on the deconvolved image, and agrees very well with Zemax simulations, as well 
as the relative intensities of the secondary lobes with respect to the central peak. On the over-resolved PSF, we note the 
presence of secondary dark rings in the Y direction, as well as pronounced diffraction spikes in the X and Y directions, 
on top of the normal JWST diffraction pattern. The dark rings are artefacts caused by the lack of information introduced 
by the errors of the hexapod motions in the Y direction (see Sect. 15.2b . The root cause of the diffraction spikes is under 




Figure 6. Illustration of the gain in resolution provided by the microscanning and deconvolution method. Left: Low-resolution (native) 
5.6 //m PSF taken in the center of the field of view Middle: Corresponding reconstructed high-resolution PSF. Right: High-resolution 
Zemax simulation, to be compared with the reconstructed, over-resolved image. Note that each individual image corresponds to a 
5" X 5" region (~ 45 x 45 pixels on the LR images, 318 x 318 pixels on the high-res images). The ghost spot below the PSF is due 
to a reflection in the window of the cryostat. All images are displayed on a logarithmic scale to exhibit the secondary lobes. The same 
dynamic range was used for the display. 
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Figure 7. Low (2 first plots on the left ) and over-resolved (2 last plots on the right) PSF profiles in the X and Y directions. For the 
high-resolution PSF, the sampling is increased by a factor of 7. The data were acquired during the FM4 campaign, in the center of the 
FoV. We compare the low and high resolution PSF to low and high resolution Zemax simulations, respectively (red dashed lines). The 
Zemax simulations shown are in the nominal, polychromatic case. 



investigation, and may be due to the coating of the murors. Those diffraction features are only visible at very high signal 
to noise ratio (typically > 2000). In addition, a ghost spot is present below the PSF, due to the reflection in the window of 
the cryostat. 

The PSF images were analysed with a set of IDL routines, complemented with the use of the IRAF software. Different 
measurements of the FWHM of the PSFs have been performed (Gauss, Airy fits, radial profiles, encircled energy, radius 
containing 65% of the total energy). 

6.1 Comparison of the results of the three test campaigns (FMl, FM2 and FM4) 

The results of the PSFs measurements in the center of the MIRIM FoV are gathered in Tables [1] and |2] Table [1] shows the 
PSF characteristics for the three test campaigns carried out at CEA. Table |2] indicates the results of Zemax simulations in 
two cases: a nominal case, and a case where mechanical tolerances are included. In the latter, we select the most extreme 
deviation from the nominal case observed in the Monte Carlo draws. The Zemax simulations are polychromatic, i.e. they 
include the width of the 5.6 filter bandpass. 



Table 1. PSF measurements for the data taken during the differents test campaigns at CEA^. 



Campaign 


Method 




Data 




Low-Resolution 


Over-resolved / = 7 




Airy 


X Y 


X 


Y 


FMl 




46.93 


47.29 


R65% 


46.79 




37.96 




^cent 


0.47 




0.51 




Airy 


X Y 


X 


Y 


FM2 


53.77 51.31 


45.73 


44.04 


^65% 


42.42 




33.10 




Ecent 


0.50 




0.56 




Airy 


X Y 


X 


Y 


FM4 


54.21 50.50 


42.26 


43.61 


1^65% 


41.90 




32.64 




Eccnt 


0.51 




0.58 



"Results are indicated for the center of the MIRIM field of view. For the HR image, the sampling is increased by a factor of 7 
with respect to a native MIRIM image, i.e. one original MIRIM pixel equals to 7 pixels in the over-resolved image. Three different 
measurements are indicated: the FWHM of the 2D Airy -Bessel function- fit {X and Y directions in iim), the radius containing 65% of 
the energy in /um, and the fraction of the total energy contained in the central lobe, normalized to 5". The statistical 1 a uncertainties on 
the fitted FWHM values are of the order of 0.04 /um. 



Table 2. PSF measurents for the Zemax simulations". 



Setup Method 



Simulations 



Zemax Low-R Zemax Hi -res 



nominal 



tolerance 



X Y X Y 
^ 48.62 50.45 42.55 42.47~ 

R65% 39.14 28.59 

Ecent 0-64 0.72 

X Y X 1^ 
~49M 51M> 4426 44.89 
R65% 41.26 30.12 ~ 

Eccnt 0.61 0.68 



Airy 



"Results are indicated for two setups of Zemax simulations: a nominal case, and a case where mechanical tolerances were applied in 
the Monte-Carlos (see text for details). 



We first note that the low-resolution PSFs are systematically wider than the over-resolved ones, which is due to the 
fact that the native 5.6 /im PSF is not Nyquist-sampled. This illustrates that the microscanning technique, and associated 
deconvolution method, are crucial to verify precisely the optical quality of the imager at this wavelength. Examples of 
PSF profiles (cut in the X and Y directions) are given in Fig. |2l The comparison of the LR and HR profiles illustrates 
the gain in resolution. The better sampling of the PSF helps in reconstructing more accurately the main peak of the PSF. 
The secondary diffraction lobes are also resolved. The observed PSF are overall in good agreement with nominal Zemax 
simulations. The cores of the PSFs are also well fitted by Bessel functions (labeled as Airy fit in the Tables), although 
we note that the positions and amplitudes of the secondary lobes do not match exactly with the Airy fits (neither with the 
simulations). 

FMl test campaign exhibits a PSF and vigneting issue. An inspection of the FMl results shows that the FWHM of the 
HR PSFs are systematically larger than the measured FWHM on the HR Zemax PSF. The discrepancy is about 
-|-8 /um. This mismatch is out of the specifications of the instrument in terms of wavefront error and encircled energy. 
This issue was also revealed independently by the measurement of the MIRIM FoV, which showed a vignetting issue. 
Motivated by these tests, a careful inspection and a new metrology analysis of the instrument was carried out at CEA 
in cooperation with IAS, Orsay. These investigations concluded to an interference between the structure of the 
Three-Mirror-Anastigmatic (TMA) objective and the M4 mirror. This design problem caused a tilt of the M4 mirror 
with respect to its nominal positioning. It was observed that the rear face of mirror M4 was not flat and was tilted 
by 0.33 degree with respect to the fixation pads. The M4 mirror has been found damaged, with some bending of its 
mounting flexible pad. This tilt resulted in a degradation of the nominal PSF properties and the FoV area. 

FM2 campaign: requirements met after correction for the tilt of the M4 mirror. The cause of the interference between 
the TMA structure and the M4 mirror was removed manually. The damaged M4 mirror was replaced by the Flight 
Spare model. A second test campaign (FM2) was launched, and we performed the same microscanning analysis 
to check the optical quality of the PSF. We note that for this campaign (as well as for the FMl), the filter wheel 
assemble was not available. It was replaced by a "cyclop tool" which allows us to mount the 5.6 ^m filter The 
results show that the correction of the M4 tilt improved the optical quality. The measured FWHM with Airy fits on 
the over-resolved PSFs are of the order of 44 — 48 /^m for both X- and Y- directions, and over the whole FoV, which 
is about 2 — 4 /im wider than the HR simulated PSF. This corresponds to FWHM of 0.19 — 0.21 arcsec. In addition, 
more than 56 % of the total encircled energy within a 5" radius is contained in the first dark Airy ring (see Sect. l6.2l l. 
This is compliant with the OBA-1004 requiremen{3. 

FM4 campaign: final verification of the optical quality with the FWA. For the FM4 test campaign, the Filter Wheel 
Assembly was available and integrated into MIRIM, thus providing PSF measurements at wavelengths 5.6 — 25.5 /im. 
The measurements show a very good agreement between the reconstructed HR PSF and the HR nominal Zemax 
simulations. Fig. [8] shows the over-resolved PSF in 9 regions of the MIRIM FoV (subarrays A to I). We found that 
the OBA-1004 requirement is met in all 9 regions (see Sect. 16.21 ). At longer wavelengths (7.7 — 25.5 /im), 57 — 68 % 
of the encircled energy is contained within the first dark ring radius, which confirms that MIRIM is compliant with 
the optical quality requirements in terms of encircled energy (see Ronayette et al. 2010"*). However, we note that 
both i?659?S and i?conJl indicate that the simulated PSF with tolerances (in the worst case) is of better quality than 
the reconstructed high-resolution PSF (Tables[T]and|2|l. We are still investigating to understand why this is the case. 
The HR technique also allows us to identify some deformations of the secondary hexagonal lobe across the FoV, as 
well as a displacement of the diffraction spike. 

6.2 Encircled energy 

This section discusses further the optical quality of the MIRIM PSF in terms of encircled energy, since the OBA require- 
ments are expressed in these terms. The encircled energy E{r) of a radially symmetric point spread function 0(r) can be 

^OBA-1004: MIRI Imager data shall include > 56 % of the Encircled Energy of the image of a point source at wavelengths longward 
of 5.6 fim within the angular diameter of the first dark ring, as determined from the PSF generated using the segmented telescope model 
in the IRD. Note: The Encircled Energy is defined relative to the energy within a radius of 5 arcsec equivalent projected onto the sky. 

*i?65% is the radius containing 65% of the total energy within a 5 arcsec radius 

^-Ecent is the fraction of the total energy contained in the central lobe, normalized to 5 arcsec. 
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Figure 8. Left: Definition of the subarrays used during cold performance tests at CEA during the FM4 campaign. Right: Reconstructed 
high-resolution PSF at 5.6 /im for each of the 9 microscans performed in imager field of view (subarrays A to I). Note that each individual 
image corresponds to a zoom on the PSF equivalent to '--^ 28 x 28 pixels on the detector (200 x 200 pixels on the high-res image). 
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Figure 9. Encircled energy of the MIRI high-resolution (HR) PSF as a function of radius in arcsec (1 pixel = 0. 1 1"). The solid lines are 
the results for the nine over-resolved PSF in different areas of the imager field of view. The dashed line shows a nominal PSF simulated 
with the Zemax software. The inflection point corresponds to the first dark ring of the Airy disk of the PSF. The percentages indicate the 
fraction of the encircled energy within the first dark Airy ring radius (situated at the radius 47.4 pm). 



defined as the integral 

E{r) = 27r / (j){r')dr' (12) 
Jo 

For the purpose of comparison to the requirement document, we normalised E{r) to 1 at r = 5 arcseconds. However, our 
microscanning data acquisition was done for a sub-array of « 7" x 7" during the FMl and FM2 campaigns to speed up the 
data acquisition. To normalize the encircled energy at a radius of 5", we used the background of a full-frame PSF image. 
For the FM4 test campaign, the sub-arrays are « 14" x 14", which allows us to normalize the encircled energy directly. 

Fig. |9] shows the encircled energy of HR reconstructed PSFs in the different regions of the MIRIM FoV, as a function 
of radius in arcsec (1 pixel = 0.11"). The dotted line shows the encircled energy of the HR Zemax PSF, computed in the 
nominal case, and in the center of the FoV. The first point of inflection in E{r), i.e the value of r for which dPE/dr^ = 0, 
corresponds to the diameter of the first dark Airy ring. We find a fraction of the encircled energy of 56.1 — 59.2 % inside 
the central lobe over the whole FoV, which is compliant with the OBA-1004 requirement. 

We note that the measured variation of encircled energy with radial distance from a point image is very sensitive to the 
flux levels at large radial distances. The condition dE/dr > for all values of r must be satisfied for the measurement to 
make sense. Here, this is the case because the background values at large radial distances are high. The telescope simulator 
is at ambient temperatures, indeed producing a high thermal background. 

7. SUMMARY AND CONCLUDING REMARKS 

In this paper we have presented the first detailed analysis of the JWST Mid-InfraRed IMager (MIRIM) Point Spread 
Function (PSF) at 5.6 /im, from three cryogenic test campaigns performed at CEA. We used a high-resolution (HR) mi- 
croscanning technique to verify the MIRIM optical quality and characterize the 5.6 /im PSF with an unprecedent accuracy. 
The microscanning consists in a fine, sub-pixel scanning of a point source on the focal plane. A deconvolution algorithm is 
used to reconstruct over-resolved PSFs (by a factor of 7 times the native resolution). We summarize here our main results: 

1 . We characterized the response of the detector, which is non-linear above 20 000 ADU. The correction for the non- 
linearity improves the quality of the MIRIM PSFs by a factor of 5 %. Note that the detector and electronics used 
for these tests were not flight model hardware. A less severe non-linear response is expected on the FM detectors. 

2. The microscanning test and associated deconvolution provide a significant improvement for the characterization of 
the PSF. This test allows to resolve the diffraction pattern of MIRIM and to measure accurately the width of the 
PSF for the flight model. The FWHM and the shape of the PSF secondary lobes is as expected from HR Zemax 
simulations. 

3. The FWHM of the HR reconstructed images taken during the FMl test campaign were 4 — 6 /im wider than HR 
Zemax simulations, and the encircled energy was out-of-specifications. The root cause was a tilt of the M4 mirror. 
After correction, we show that the observed HR PSF is in excellent agreement with simulations. The FWHM is 0. 19" 
- 0.20", and 56.1 — 59.2% of the total encircled energy within a 5" radius is inside the central lobe, over the whole 
field of view, which is within specifications. At longer wavelengths (7.7 — 25.5 /im), 57 — 68 % of the encircled 
energy is contained within the first dark ring radius, which confirms that MIRIM is compliant with the optical quality 
requirements. 

4. Slight differences between the simulated PSF patterns and the MIRIM HR reconstructed data are exhibited. Some of 
these features (dark arcs) are associated with the limits of the deconvolution method itself, in particular due to errors 
in the relative positions of the images and the associated lack of information. The other features (deformation of the 
secondary lobes and diffraction spikes) require further investigation to understand their root cause. 

The microscanning technique and deconvolution algorithm described in this paper have been particularly powerful 
to verify the optical performance of MIRIM and identify defects in the instrument. This high-precision characterization 
of the PSF at the shortest operating wavelength will be of importance for the photometric calibration of MIRJ, but also 
for possible scientific applications, like point source extraction or map-making, where a precise knowledge of the PSF is 
required. This technique may also be applied for the MIRI Medium Resolution Spectrometer (MRS) integral field unit. In 



addition, we note that the dithering mirror implemented on the JWST may be used to perform over-resolution, following 
the methods described here to deconvolve the images. This may push even further the capabihties of the JWST in terms of 
angular resolution. 
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